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Abstract 

A rarefied-flow shuttle aerodynamic model spanning the 
hypersonic continuum to the free molecule-flow regime has been 
formulated. The model development has evolved from the High 
Resolution Accelerometer Package (HiRAP) experiment conducted on 
the Orbiter since 1983. The complete model is described in detail. The 
model includes normal and axial hypersonic continuum coefficient 
equations as functions of angle-of-attack, body flap deflection, and 
elevon deflection. Normal and axial free molecule flow coefficient 
equations as a function of angle-of-attack are presented, along with 
flight derived rarefied-flow transition bridging formulae. Comparisons 
are made with data from the Operational Aerodynamic Design Data 
Book (OADDB), applicable wind-tunnel data, and recent flight data 
from STS-35 and STS-40. The flight-derived model aerodynamic force 
coefficient ratio is in good agreement with the wind-tunnel data and 
predicts the flight measured force coefficient ratios on STS-35 and STS- 
40. The model is not, however, in good agreement with the OADDB. But, 
the current OADDB does not predict the flight data force coefficient 
ratios of either STS-35 or STS-40 as accurately as the flight-derived 
model. Also, the OADDB differs with the wind-tunnel force coefficient 
ratio data. 


Nomenclature 

a =acceleration measurements 

C =aerodynamic force coefficient 
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c 

=normalized aerodynamic coefficient 

Kn 

=Knudsen number: mean-free path/12.058 m 

M 

=Mach number 

V 

=velocity 

Vbar 

=viscous parameter 

5 

=deflection angle 

a 

=angle-of-attack 

Acronyms 

HiRAP 

=High Resolution Accelerometer Package 

IMU 

=Inertial Measurement Unit 

OADDB 

=Operational Aerodynamic Design Data Book 

STS 

=Space Transportation System 

Subscripts 

A 

= axial force 

bf 

=body flap 

c 

=continuum 

EL 

=elevon 

f 

=free molecular flow 

N 

^normal force 

a 

=angle-of-attack 

x,y,z 

=body axes 


Introduction 

Prediction of vehicle aerodynamic coefficients come from several 
sources, e.g., wind-tunnel testing, computational simulations, and 
aerodynamic models. Historically, starting with the Wright brothers, 
when appropriate scaling parameters) are identified, wind-tunnel data 
have successfully been applied to predicting aerodynamic flight 
behavior. But, with the advent of space flight, and the interest in 
hypersonic rarefied-flow phenomena, current tunnel capabilities are 
being pushed to the technology limits for providing flight predictions. 
Over the last decade or so, a parallel effort for providing wind-tunnel 
like information has been developed using computer simulations of the 
flow-field. These computational methods are becoming more potent 
with continued computer innovations and the theoretical bounds appear 
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limitless. However, there are technical problems which are more 
economically served by a simpler approach, for instance, preflight 
design studies. In these cases, a simple aerodynamic model is adequate. 
This report provides an aerodynamic force model for applications on 
winged entry vehicles in the hypersonic rarefied-flow flight regime. 

In the recent past, the Shuttle was especially instrumented with a 
variety of experiments to measure aerodynamic and 
aerothermodynamic parameters. 1 This capability provided an 
opportunity to make repeated in-situ measurements in the rarefied- 
flow regime so that the predictive models could be developed with flight 
data. A rarefied-flow aerodynamic model has been formulated from the 
High Resolution Accelerometer Package (HiRAP) experiment.2 By 
examining the force coefficient ratio on multiple flights, the rarefied- 
flow transition aerodynamic functional behavior has been extracted. 
This provided the keystone to a simple model which can be used to 
predict rarefied-flow aerodynamics of a winged reentry vehicle, e.g., the 
Shuttle. The model incorporates aerodynamic coefficient derivative 
data from the Operational Aerodynamic Design Data Book ( OADDB )3 
at the rarefied-flow transition boundaries, that is in the hypersonic 
continuum and in the free molecule-flow regime. In the hypersonic 
continuum regime, the absolute value of the individual force coefficients 
in the model have been adjusted to values determined on STS-61C, 
when hypersonic continuum pressure measurements first became 
available. The Orb iter control surface aerodynamic contributions, i.e., 
the body-flap and elevon coefficients, are generated from curve fits to 
the OADDB data. At the other end of the rarefied-flow transition 
boundary, the free-molecule flow force coefficient values are also 
generated from curve fits to the OADDB data. 

The purpose of this report is to provide a complete description of the 
flight model and to show its usefulness to predict the aerodynamic force 
coefficients of the Orbiter in the rarefied-flow regime. To demonstrate 
this, the model is used to predict the force coefficient ratios of two recent 
flights, STS-35 and STS-40. For these flights, a similar comparison is 
made using the OADDB. Also, both the flight model and the OADDB are 
compared with rarefied-flow wind-tunnel data taken early in the 
Shuttle vehicle development. 
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Model Description 

The sign convention used for the model is the standard used in 
the Orbiter program documentation and is presented in Fig. 1. 

The OADDB data come from the "change 3" update which reflect 
the most recent normal and axial coefficient Shuttle flight information. 
A subset of OADDB data were chosen which includes parameters 
appropriate to the comparisons discussed later. These parameters 

include: angle of attack with range of 0° to 60° for the payload bay doors 
closed (corresponding to free molecule-flow conditions), angle of attack 

with range of 35° to 45° for the entry region, elevon deflection with 

range of -15° to 22.5°, body flap deflection with range of -11.7° to 22.5°; 
viscous parameter with range of 0.001 to 0.05, and Mach number with 
range of 15 to 30. The entry coefficients were obtained by summing 
contributions from the basic, high altitude, real gas, viscous interaction, 
body flap, elevator/aileron (elevon), and vehicle components. 

The rarefied-flow Shuttle aerodynamics flight model discussed in 
this report is a simple model. It contains three segments: the 
hypersonic continuum, the free molecule flow, and the transition 
bridging formulae. Each segment is described in detail. 

Hypersonic Continuum 

The formulae for the hypersonic continuum normal and axial 
coefficients as a function of a are: 


C =-925704 X10 5 oc 2 +523808 X10 2 a-£39782 
N 

c, a 

C =5 £6689 X10 -7 a 3 -6.72027 Xl0 _5 a 2 +3 32044 XlO" 3 a- .0086314 
A 

c/u. 

These equations are for the Shuttle vehicle at an angle-of-attack 

envelope of 35° < a < 45°. The above coefficient functions are resulting 
curve fits to the basic data in the OADDB with adjustments resulting 
from an earlier analysis using flight data from mission 61C. 2 Unlike the 
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more sophisticated OADDB, the flight model basic coefficients previously 
given do not separate physical phenomena such as, real gas and viscous 
effects . In effect, each equation results in one combined coefficient 
allowing for computation ease and, as discussed later, not introducing 
significant errors. Figure 2 shows the calculations using the above 

equations for angle-of-attack between 35° and 45°. Also included on the 
graphs are a comparison of the flight model force coefficients angle-of- 
attack behavior with the OADDB data. The conditions for the OADDB 

data are M=15, Vbar=0.005, and &bf=5 EL =0°. These conditions 

corresponds to an altitude of approximately 58 km. Two OADDB data 
sets are shown; one is the OADDB basic coefficients (the circles), while 
the other (the squares) has the real gas, viscous and vehicle effects 
included, i.e., a composite coefficient, comparable to the above equations. 
The primary difference between the individual OADDB basic 
coefficients and the points labeled, "OADDB (total)" is real gas effects for 
C N and viscous effects for C A . Comparing these OADDB composite 

coefficients with the flight model, it is seen that the flight model C N 

coefficient is consistently lower by about 7 percent than the OADDB, 
while the C A coefficient is slightly higher (about 2 percent) for all 

angles-of-attack. These differences will produce noticeable differences 
(about 9 percent) in the coefficient ratio, to be discused later. 

The Shuttle control surfaces modeled are the hypersonic 

continuum body flap and elevons at a = 40° and M = 20. The 
normal and axial contributions for each are: 

BODY FLAP 


for -11.7° < 5. < 22.5°, 
bf 

C r = 4.46278 x\0~ 4 + 1.92931 xlO _3 5, + 3.6029 xl(T 5 8?, 

N bf bf 

c,bf 

C A = 2.39178 x 10 -4 + 4.81862 x lO^S, , + 3.03937 x 10“ 5 S^ - 3.500869 x 10 -7 8? 
A L , bf bf k 
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ELEVONS 


for -15° < 8 < 15 u , 

EL 


s 0 


C = 1.00333 xlO -3 +5.76021 x 10 3 S r ,+ 1.18538 x 10 4 5 2 
N EL EL 

c,EL 

C = -1.3981 xlO _3 + 1.55864 xlO - ^, + 8.6081 xl0 _5 5 2 
A EL EL 

c,EL 

Figures 3 and 4 show the above functions compared with the data 
from the OADDB for the body flap and elevons, respectively. 

Then, the total hypersonic continuum coefficients are: 


C N ~ C N +C N +C N 
c c, a c,bf c,EL 

C A =C A +C A +C A 

A A A A 

c c, a c,bf c,EL 


Free Molecule Flow 


The formulae for the free molecule flow normal and axial 
coefficients as a function of a are: 


C = -7.16528 xlO _6 a 3 + 9.66197 x lO^a 2 + 9.18422x10 3 a + 1.58739x10 3 
N 

/.a 

C = -1.171 17 x 10 _5 a 3 + 5.92205 x 10 _4 a 2 +.0164864a+.751 105 

A f , « 


The above free molecule flow normal and axial coefficient 

equations are for the Shuttle (payload doors closed) at 0° < a < 60°. 
A comparison of these functions with the data from the OADDB is 

shown in Fig. 5. The body flap and elevons contributions, that is, 


'N 


N 


f,bf f,bf f,EL f,EL 
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in the free molecule flow regime are undetermined and are 
currently set to zero. 

Then, the total free molecule flow coefficients are: 


r — c + c 

N N N ' ~N 

f f.a f,bf f,EL 


+ C. 

+ c 


C. = C A +C A 

A A A A 

f f, a f,bf f,EL 


Bridging Formulae 

The method to bridge between the hypersonic continuum to the free 
molecule flow regime includes the following "bridging" formulae^ 


C N = exp| 


-0.2998^ 1.3849- log, n Kn 
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1.7128 


otherwise C, 7 =1.0 

N 


C A =exp 
A 


-0.2262^1.2042 - log i{) Kn^j 


1.8410' 


if log iQ A>i < 1.3849 


if IoQjq Kn < 1.2042 


otherwise C, =1.0 

A 


where the normalized normal and axial coefficients use Knudsen 
number, Kn, as the independent parameter. The normalized 
aerodynamic coefficients are used to generate the actual cofficients by 
the following: 






c 







c 


Flight Data 


The flight data for STS-35 and STS-40 is summarized in this 
section. The Orbiter position, velocity, and orientation information 
is obtained from a "best estimate trajectory" process.^?® This 
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process includes using flight measurements from the IMU and is 
the source of the IMU acceleration ratios shown later. Figure 6 shows 
the Orbiter velocity and corresponding Kn as a function of altitude for 
STS-35 and STS-40. The Kn is computed using the 1976 U.S. standard 
atmosphere7 and a reference length of 12.058 m, the Orbiter's mean 
aerodynamic chord. The Kn parameter serves as a guide for the vehicle 
flow regime. Namely, large Kn (roughly > 10 ) indicate free molecule 
flow behavior with few intermolecular interactions, while small Kn 

(roughly < 10' 3 ) indicate near continuum flow where intermolecular 
collisions are an important consideration to the flow physics. Clearly, 
most of the altitudes included in Fig. 6 are for the rarefied-flow 
transition regime, where the vehicle environment is "transitioning" 
from free molecule flow to continuum flow. In addition to a flow 
indicator, Kn is also the independent parameter for the bridging 
formulae presented earlier. 

Figure 7 gives the angle-of-attack profile of the Orbiter vehicle for 
STS-35 and STS-40. Throughout the entire rarefied-flow transition into 
the hypersonic continuum region, about 170 km to 60 km, the angle-of- 

attack is close to 40°. The corresponding control surface settings for 
these flights are shown in Fig. 8 and Fig. 9. Both body flap and elevon are 
held constant in their stored positions until shortly below 90 km, when 
they are activated, mostly in the pitch direction, in order to establish 
longitudinal control. 

An application of the STS-35 and STS-40 flight data with the 
rarefied-flow aerodynamics flight model described in the preceding 
section produces the results shown in Fig. 10. These results show 
angle-of-attack effects which influence both the and Cj^ coefficients. 

The angle-of-attack variations are easily seen in C N at both low and 
high altitudes (see Fig. 7 for comparisons). Due to scale selection, the 
angle-of-attack variations are not easily seen at low altitudes, but are 
discemable at high altitudes. 


Model Comparisons 


With Flight Data 

A method of examining aerodynamic models with flight data is to 
compare the model coefficient ratio with the corresponding 
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acceleration measurement ratio since, 

Cn _ az 
C A a x 

The use of these ratios alleviates the need to determine dynamic 
pressure and thus, density which is a widely varying quantity. 2 
Figure 11 shows 2 graphs, both of which contain measurements 
of the ratio of acceleration from the IMU and HiRAP on STS-35 
(the IMU data are below 97 km). Figure 11(a) shows the results 
of application of the aerodynamic coefficients using the OADDB. 

Clearly, there are differences of the coefficient ratio throughout the 
rarefied-flow regime with the most notable deviations, about 10 percent, 
in the hypersonic continuum. Figure 11(b) includes a comparison with 
the aerodynamics generated with the flight model and is labeled "flight 
model." The difference between this model and the flight data is very 
small, less than 1 percent throughout the entire regime. Fig. 12 is a 
similar set of comparisons using data from STS-40. Figure 12(a) shows 
again the comparison of the flight data with the OADDB, while Fig. 

12(b) shows a comparison with the flight model. The results further 
confirm STS-35 comparisons, namely, the OADDB continuum 
aerodynamics do not satisfactorily match the flight data. 

With Wind-Tunnel Data 

Wind-tunnel data were collected by the Orbiter Project on a 1.0 
percent scale model in a hypersonic shock tunnel facility^ in order to 
examine viscous interaction effects. The test conditions were Mach 
numbers up to 16 and viscous interaction parameters values up to 0.06. 
Shown in Fig. 13 are the C N and C A test results, in terms of a ratio, for 

a= 40° and 8^ = 5 EL = 0°, as a function of Kn. Also shown in the figure 
are the OADDB data and the flight model. As seen, the ratio of wind- 
tunnel force coefficient data is in excellent agreement with the flight 
model. This is somewhat fortunate, since clearly combinations of 
viscous and real gas effects must cancel in the ratio for these coefficients. 

Also shown in the figure are the ratio data from OADDB which show 

-3 

differences with the wind-tunnel data. (Note, at Kn=10" , corresponding 
to about 86.868 km (300,000 ft), there is a small discontinuity in the 
OADDB.) The differences between the OADDB data and the wind- 
tunnel data are small, but appear statistically significant insofar as the 
functional behavior of the transition into the hypersonic continuum flow 
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is not being followed. This is not true for the flight model which shows 
good agreement with the tunnel data. This fact, coupled with the good 
agreement with the STS-35 and STS-40 flight data, provides an 
argument that the flight model has good predictive capabilities 
throughout the entire rarefied-flow regimes. Further, its application is 
relatively simple. 


Summary 

A rarefied-flow Shuttle aerodynamics model with data from the 
OADDB and HiRAP flight measurements has been formulated. This 
model is compared with flight data from STS-35 and STS-40, wind 
tunnel data, and the current OADDB. The method of comparison uses 
the force coefficient ratio in order to eliminate the dependency on 
density. There is excellent agreement between the flight aerodynamic 
model and the flight data( ~1 percent) and the wind-tunnel data when 
expressed as a ratio. Similar comparisons with the OADDB exhibit a 
problem as the hypersonic continuum is approached. Namely, the 
OADDB hypersonic continuum coefficient ratio values show a 
consistent significant difference with the flight data (about 10 percent) 
and the wind-tunnel data. 
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(b) for aerodynamic control surfaces. 


Fig. 1 Sign conventions. 
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Fig. 2 Body axes coefficients comparison - hypersonic continuum. 
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Fig. 5 Body axes coefficient comparison - free molecular flow. 



Fig. 6 Orbiter flight conditions. 
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Fig. 7 Orbiter angle-of-attack flight profiles. 



Fig. 8 Orbiter control surface settings for STS-35 
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Fig. 9 Orbiter control surface settings for STS-40. 
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Fig. 1 1 Comparison of predictions with STS-35 measured acceleration 
ratios. 







Fig. 12 Comparison of predictions with STS-40 measured acceleration 
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